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ABSTRACT 

Individual gamma ray bursts (GRBs) have very diverse time behavior - from a single 
pulse to a long complex sequence of chaotic pulses of different timescales. I studied 
light curves of GRBs using data from the Compton Gamma Ray Observatory's Burst 
and Transient Source Experiment (BATSE) and found that the average post-peak time 
history for a sample of 460 bursts obeys an unique and simple analytical law: I ~ 
exp( — (t/to) 1 / 3 ) where t is time measured from the peak of the event and t is a constant 
ranging from 0.3 sec for strong bursts to ~ I sec for weak bursts. The average peak 
aligned profile follows this law with good accuracy in the whole range availible for analysis 
(from fractions of a second to ~ 150 seconds after the peak). Such a law with a single 
time constant characterising the overall sample of GRBs should have important physical 
meaning. The same analysis for solar flares gave a dependence close to exp( — (t/to) 1 ^ 2 ) 
and some indication of a tendency to change the power index to 1/3 at t > 50 sec after 
the peak of the flare. The dependence of to versus brightness of GRBs is presented. 
The fact that to depends on the brightness apparantly confirms the recently discovered 
effect of time dilation of weak bursts which has a possible cosmological interpretation. 
The time dilation is detected at a confidence level of 7a and it is slightly larger than 
was previously reported. The statistics (460 GRBs), however, still does not allow us 
to rule out an intrinsic anticorrelation between intensity and duration of GRBs as a 
possible explanation of the effect. Probably it is possible to confirm (or disprove) the 
cosmological origin of GRBs using the total existing statistics of strong bursts. 

Subject headings: cosmology: theory - gamma rays: bursts 
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1. Introduction 

The characteristic distance to sources of gamma 
ray bursts is still unknown. A possibility that GBRs 
could originate from cosmological distances stimu- 
lated studies of their light curves in order to find 
signatures of the time dilation of weak GRBs rel- 
ative to strong GRBs associated with a difference 
in their respective redshifts (Paczyhski 1992; Piran 
1992). These studies used different techniques. Davis 
et al. (1994) measured widths of individual pulses in 
multipeaked GRBs. They found a time dilation of 
order 2 for pulses in weak GRBs. Mitrofanov et al. 
(1994) compared average profiles of weak and strong 
GRBs and did not find any significant time dilation in 
a sample of 260 GRBs from the first BATSE catalog. 
Norris et al. (1994) using the same method as Mitro- 
fanov et al. (1994) as well as other tests found a time 
dilation by a factor 2 using the larger statistics of the 
second BATSE catalog. Finally Norris et al. (1995) 
confirmed the time dilation effect measuring duration 
distributions of bright and dim GRBs. 

All these studies did not consider possible charac- 
teristic shapes of measured distributions. One can 
make better estimates of the time dilation if one fits 
the data with a meaningful parametrized expression 
for the shape of distributions. In order to study pos- 
sible signatures of cosmological effects in the time be- 
haviour of GRBs, I decided to investigate the aver- 
age shape of their light curves. To define the average 
time history one needs a time reference point for each 
GRB, the most reasonable being the time of the peak 
of the event. With this reference time one obtains an 
average peak aligned profile of GRBs, (Mitrofanov et 
al. 1994; Norris et al. 1994). I adopted the same 
approach and concentrated on the post-peak slope of 
the average profile as the pre-peak slope is influenced 
by the triggering effects and the corresponding data 
have a worse quality (the pre-peak slope is steeper, 
see Nemiroff et al. 1994). The purpose was to find an 
analytical expression that fitted the profile. 

2. The average peak aligned profile of GRBs 

All data utilised in this work were extracted from 
the publically available BATSE database in God- 
dard Space Flight Center. Time profiles of GBRs 
are derived from 64 ms resolution Large Area De- 
tectors data by summing up counts in 4 discrimina- 
tor channels covering the 25 keV - 1000 keV energy 
band. For relatively bright GRBs (peak count rates 



P > 4000s -1 ), positions of the peak were determined 
by the brightest 64 ms bin. For weaker events the 
peak was determined to be at the center of the 128, 
256 or 512 ms brightest interval (in order of decreas- 
ing P). 

To subtract the background I used linear fits - this 
procedure is not only simple but is also optimal be- 
cause it preserves the statistics of weak tails of GBRs 
better than quadratic or higher order polynomial fits. 
To prepare the fitting procedure I made a visual ex- 
amination of all GBR's time profiles determing fitting 
windows in a way to avoid regions with probable GBR 
contributions. For many events the 1024 ms resolu- 
tion data covering a 470 s time interval were also used 
to determine two fitting windows - one before and one 
after a burst. Typically, the right fitting window was 
separated by at least 150 s from the burst peak (ex- 
cept for apparently short bursts). 

After discarding all events with severe data gaps 
and a couple of events with a patalogically variable 
background I found 460 useful GBRs in the second 
BATSE catalog. 

To divide GBRs into brightness groups I used 64 
ms and 1024 ms integration time criteria for the peak 
brightness. This work is mainly based on the tradi- 
tional 1024 ms brightness selection timescale because 
of its compromise between peak brightness and total 
energy release criteria. I did not select bursts using 
their durations as a criteria. The resulting profiles 
therefore include all GRBs, even those which are con- 
fined to within a single 64 ms bin. This is impor- 
tant because removing short GRBs from the sample 
would deform the average profile near its peak. Un- 
resolved events give the proper normalization for the 
zero-point of the profile. 

The average peak aligned time profiles obtained 
with this method extend for more than 150 seconds 
after the peak. I found that the shape of the average 
profile obeys with good precision a simple law: 

I(t)/P = exp(-(*Ao) 1/3 ) 

where t is the time after the peak, P is the peak count 
rate (intensity) and to is a constant. This is demon- 
strated in Figure 1 where the average time profiles 
for two brightness groups are plotted as log(J) vs t 1 ! 3 . 
The time constant to depends on the brightness group 
(to = 0.48 s for the bright group, to = 1.01 s for the 
weak group) apparantly demonstrating time dilation 
of weak bursts in approximate agreement with results 
of Norris et al. (1994, 1995). 
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The time profile averaged over 409 GBRs (the 
weakest events were excluded from this sample) has 
an intermediate to (to = 0.72 s) and closely follows the 
exp(— (t/to) 1 / 3 law. However, if we believe in the cos- 
mological nature of the time dilation it makes sense 
to rescale all light curves to the same redshift. This 
procedure allows us to get rid of a possible slight cur- 
vature resulting from mixing linear profiles of differ- 
ent slopes. I calculated the time profile averaged over 
409 GRBs, where time profiles of 3 brightness groups 
covering the whole sample were rescaled to to = 1- 
The residuals of the linear fit to this time profile are 
presented in Fig. lb, while the time profile itself is 
presented in Fig. 2. The linear behavior is clearly 
visible up to t ~ 200 s. The deviation from the linear 
fit does not exceed 12% within the first 100 s, over 
which the profile amplitude varies by two orders of 
magnitude. 

The result of the two parameter linear fit, If = a + 
bt, in the interval 0.5 < t 1 ! 3 < 5 is: b= 1.07, a = 1.35 
(the standard \ 2 errors are meaninglessly small in this 
kind of fit due to correlations between bins), with a 
corresponding to the value of the profile extrapolated 
to t = 0. The extrapolated value exceeds the actual 
profile value at t = which is 1 by definition. This 
difference is mainly caused by the finite width of the 
first bin (see Fig. lb) 

In order to test the exp(— (t/to) 1 ^ 3 ) law, I varied 
the power index a in a exp(— (t/to) a ) law, rescaled 
the abscissa to t a , and repeated the linear fit (the 
profile was represented by only 10 bins in this test 
in order to avoid a strong correlations between bins). 
The result of a X 2 test is a = 0.344±0.025 at the 90% 
confidence level, which is sufficiently close to 1/3 that 
it is reasonable to suppose that a could be exactly 
1/3. 

I checked the behavior of the peak aligned pro- 
file for solar flares using an arbitrary sample of 290 
flares from the same BATSE database. Fig. 2 
presents the average time profile of solar flares to- 
gether with the general profile of GRBs in coordi- 
nates: — log(— log(7/7o)) vs log(t), where 7o is the 
average intensity of the profile extrapolated to t = 
(it slightly differs from I p as mentioned above). 
These coordinates are convenient to demonstrate de- 
pendences such as exp(— x a ) in the same way as a 
log — log plot demonstrates power laws x a . The aver- 
age profile of solar flares behaves as exp(— (t/to) 1 ^ 2 ), 
where to ~ 6 sec for the first 40 seconds after peak, 
then it changes power index to the same value, a = 



1/3, as for GRBs. Unfortunately, the data cover a too 
short time interval to make this a strong conclusion 
(I used only 64 ms resolution data covering 240 sec- 
onds, 1.024 sec resolution data may extend the curve 
somewhat). Independently of whether the profile of 
solar flares actually change index to 1/3 or not, peak 
aligned time profiles of GRBs and of solar flares have 
the same type of dependences. 

If one looks at time profiles of individual GRBs 
plotted as log(7) vs t, one finds well shaped expo- 
nential tails of the pulses with the decay time con- 
stant varying from milliseconds to approximately ten 
seconds. However, the exp(— (t/to) 1 ^ 3 ) behavior is 
hardly a sum of these exponentials. Rather it seems 
that the main contribution to the average time pro- 
file, decaying slower than an exponential, comes from 
multiple pulses within individual GRBs. 

I believe that the observed shape of the time profile 
could be important for understanding GRBs. How- 
ever even before the exp(— (t/to) 1 ^ 3 ) law is under- 
stood, the knowledge of average GBRs time behavior 
presents an opportunity to measure the effect of time 
dilation of weak GRBs with better precision than was 
possible before. 

3. Time dilation 

The exp(— (t/to) 1 ^ 3 ) shape of the time profile of 
GRBs is a property which complicates the measure- 
ment of the time dilation effect. When fitting the 
profile with a straight line in log(7/7 p ) - t 1 ! 3 coordi- 
nates, the resulting slope is a cubic root of the time 
constant and all errors are multiplied by 3 when con- 
verting to normal time scale. This problem cannot 
be bypassed by rescaling to linear time and it also 
complicates the analysis of distributions of durations 
(the dispersion of durations is too large). This loss of 
accuracy is a kind of invariant and this is the proba- 
ble reason why Mitrofanov et al. (1994) did not find 
any time dilation effect. However, knowing the ac- 
tual shape of average profiles we can not only detect 
the time dilation but also to measure the dependence 
of the time dilation on the brightness of GRBs with 
satisfactory accuracy. 

I measured to for different brightness intervals us- 
ing linear \ 2 fits in log(7/7 p ) - t 1 ! 3 space. The results 
are presented in Fig. 3. It includes two sets of sam- 
ples - one set with wide brightness intervals (7 sam- 
ples, 70 events in each sample on average, except for 
the brightest sample which includes only 31 events) 
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and a second set with narrow brightness intervals (23 
samples, 20 events in each on average). The pur- 
pose of the narrow brightness range sampling was to 
get an independent estimate of errors. The errors 
estimated from the \ 2 fiting procedure are not suffi- 
ciently reliable (usually being underestimated) due to 
correlations between neighboring time bins. The rms 
errors were extracted from root mean square devia- 
tion of to iteration power law fit to 23 narrow bright- 
ness range samples. These errors were then applied 
to wide brightness range samples by reducing the er- 
rors assuming that the relative rms error of a sample 
is proportional to TV -1 / 2 , where N is the number of 
events in the sample. 

The time dilation effect is slightly larger than found 
by Norris et al. (1994, 1995). Relative to the brightest 
sample, the time dilation i.e. the average ratio of to 
for the brightest sample to to for 3 weakest samples, 
is 3.2 ± 0.7 (The error is associated mainly with the 
brightest sample which includes only 31 events). 

Fig. 3 indicates that the dependence of to might 
be flatter in the dimmer part of the peak inten- 
sity range than in the brighter part (in the cosmo- 
logical interpretation this could mean that all weak 
bursts are near the edge of their redshift distribu- 
tion). However the statistical significance of such a 
flattening is small and the distribution can be fit- 
ted with a simple power law with a satisfactory \ 2 ■ 
I performed a power law fit to the wide brightness 
range samples using the rms errors obtained in the 
procedure described above. The result of such a fit is 
[3 = —0.24 ± 0.035 where [3 is the power index in the 
dependence t = (0.36 ± 0.04)s x (7f 024 /60000s- 1 )-' 3 , 
7f 024 is the peak intensity on the 1024 ms integration 
time scale in units of counts per second. Note that the 
time dilation effect is detected in terms of a nonzero 
[3 at a 7cr level. 

These results were obtained using count rates for 
a wide energy band 25 - 1000 keV. However, GRBs 
are known to have shorter durations at high energies 
relative to low energies (e.g., Mitrofanov et al. 1994). 
I checked the dependence on energy band measuring 
to for BATSE channels 1 and 2 (25 keV - 100 keV) and 
channels 3 and 4 (100 keV -lOOOkeV) separately. The 
dependence is moderate: the to for the softer energy 
band is 1.5 times larger than for the hard energy band 
for both bright and dim bursts. 

The dependence of to on energy band means that 
if we adopt the cosmological interpretation, the mea- 
sured time dilation effect is slightly underestimated. 



Spectra of weak bursts should be redshifted (and they 
actually are softer, Paciesas et al. 1991, Mitrofanov 
et al. 1994) and a soft longer component of bursts 
would move below threshold. Thus what we mea- 
sure for dimmer bursts is their shorter component. 
The possible correction for this effect is +25 — 1-40 % 
(taking into account the width of the energy band of 
channels 1 and 2). The index [3 then becomes close to 
0.3 and the maximal measured time dilation should 
be corrected to ~ 4 ± 0.9. Of course, this correction 
makes sense only in the cosmological scenario. 

4. Discussion 

The average peak aligned profile of GRBs obeys a 
exp(— (t/to) 1 / 3 ) law for a wide range of durations and 
amplitudes with good accuracy. This could be due to 
something more fundamental than a coincidence re- 
sulting from adding exponentials with different decay 
times. As I already mentioned, the main contribution 
to the average time profile comes from chaotic repeat- 
ing pulses rather than the exponential decays of main 
pulses. This means that the law should be an intrin- 
sic feature of some stochastic process responsible for 
the time behavior of GRBs. One possible approach to 
the problem is to study a number of simple stochastic 
toy models in the style of Rosner & Vaiana (1978). 

It is also reasonable to suppose that the similarity 
in time behaviour between solar flares and GRBs is 
not a coinsidence. It was pointed out in many works 
(e.g. Dennis 1988, Hurley 1989) that solar flares and 
GRBs have similar temporal and spectral properties. 
One can conclude from Fig. 2 that GRBs just is a 
cleaner variant of the same class of behavior. This 
could mean that GRBs are associated with the same 
class of processes as solar flares i.e. magnetic recon- 
nection. Magnetic reconnection can in principle be 
responsible for GRB's light curves not only in the case 
of a galactic origin of GRBs. A cosmological scenario 
necessarily includes a relativistic expanding fireball of 
a relativistic jet (e.g. Piran 1994). Both the fireball 
and the jet can easily be magnetically dominated, and 
an essential fraction of the total energy can be carried 
by the magnetic field. If the field is turbulent (there 
exist models of fireballs with turbulent field growth, 
see Meszaros & Rees 1994) it will reconnect during 
the expansion releasing energy and the light curve of 
a GRB will just reflect the history of these reconnec- 
tions. 

Concerning time dilations I believe that at the 
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present statistical level they are consistent both with 
the cosmological scenario and with the existence of an 
anticorrelation between intensity and duration. The 
only way to confirm or to rule out the cosmological 
interpretation of time dilations is to use a larger sta- 
tistical sample. The important point is that probably 
the necessary amount of statistics already exists in 
old archives. 

The cosmological scenario implies a z ~ region 
where there are no time dilations (unless there existed 
a strong evolutionary effect that completly suppressed 
GBRs in this region). This means that the time di- 
lation curve should be flat for the brightest range of 
GRBs. The estimated time dilation is already large 
and probably the brightest sample in Fig. 3 should 
be close to z ~ region. 

The statistics which I utilized is too poor to con- 
firm this. It covers only 2 years of observations while 
the total observation time of GRBs with all instru- 
ments (PVO, WATCH-GRANAT, KONUS etc.) is a 
few times larger. All data since the late 1960s could 
be utilized because the sensitivity of the instrument 
does not matter in this case as only light curves of 
strongest bursts are of interest. If a reliable signa- 
ture of a z ~ region will be detected in the depen- 
dence of to vs peak intensity, it would leave almost no 
doubt that the cosmological interpretation is correct. 
If, on the contrary, the strongest bursts demonstrate 
a monotonic decrease of the to then the time dila- 
tion between strongest and weakest bursts will exceed 
some reasonable cosmological limit and the effect of 
time dilations must require some other explanation. 
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Fig. 1. — (a) Average peak aligned profiles for two 
brightness groups (/f 24 > 7800s -1 , 95 GRBs, thick 
solid curve, and 860s" 1 < J[ 024 < 3000s" 1 , 218 
GRBs, thin curve), (b) Residuals of the linear fit to 
the average profile for 409 GRBs rescaled to to = 1- 
Dotted line indicates the la envelope. The straight 
line fragment on the left represents the first 64 ms bin 
assuming the count rate distribution within the bin 
to be uniform in t 1 ! 3 . 
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Fig. 2. — Average peak aligned profiles for 409 GRBs 
and for 290 solar flares in coordinates chosen to dis- 
play an exp(— t a ) dependence as a straight line with 
slope a. Jo is the profile intensity extrapolated to 
t = 0. This removes uncertainties associated with 
the finite width of the first bin. The cutoff at large t 
corresponds to the 240 s duration of data records. 
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Fig. 3. — Time constant, to versus peak intensity, 
Jf 024 , in log - log coordinates. Plus signs represents 
to for samples of wide brightness ranges. The rms 
errors were estimated from the dispersion of to for 
samples with narrow brightness ranges [thin lines). 
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